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We study carrier transport through graphene on SrTiC>3 substrates by considering relative con- 
tributions of Coulomb and resonant impurity scattering to graphene resistivity. We establish that 
charged impurity scattering must dominate graphene transport as the charge neutrality point is 
approached by lowering the carrier density, and in the higher density regime away from the neu- 
trality point a dual model including both charged impurities and resonant defects gives an excel- 
lent description of graphene transport on SrTiOa substrates. We further establish that the non- 
universal high-density behavior of o~(n) in different graphene samples on various substrates arises 
from the competition among different scattering mechanisms, and it is in principle entirely possible 
for graphene transport to be dominated by qualitatively different scattering mechanisms at high 
and low carrier densities. 



I. INTRODUCTION 

The substrate material has significant influence on 
graphene transport properties, which has been confirmed 
by many different experimental groups^—. The exten- 
sively used model^i of Coulomb disorder in the environ- 
ment has earlier been found to provide a reasonable the- 
oretical description for graphene transport on Si02 sub- 
strates from many different groups^ as well as on several 
other substrates^^. In a recent paper—, however, Couto 
et al. claim that the long-range Coulomb impurities do 
not play an important role based on their measured con- 
ductivity, cr(n), of graphene at high carrier densities on 
SrTiOa substrates. In addition, they argued that their 
transport data can be explained quantitatively by a so- 
called "resonant scattering" modelii, which gives the fol- 
lowing expression for the carrier density (n) dependence 
of the conductivity a: a{n) = 2%. ^- In 2 (V nnR 2 ) , where 
Hi, R are respectively the concentration and the range of 
the resonant scattering defects in graphene. Motivated 
by this experiment^, we revisit the question of the role 
of various types of disorder on different substrates in con- 
trolling the density-dependent conductivity of graphene, 
comparing, in particular, the low- and high-density con- 
ductivity limited by distinct scattering mechanisms. In 
particular, the low-density graphene transport behavior 
on SrTiC>3 (or for that matter, on any other substrate) 
is likely to be always dominated by long-range Coulomb 
disorder, independent of the other scattering mechanisms 
(e.g. short-range disorder, resonant scattering) which 
may be operational at high carrier density. One impor- 
tant finding of our current work is to show that it is 
indeed possible for graphene transport to be determined 
by different scattering mechanisms at low and high car- 
rier densities, a result which has been implicit in earlier 
works in this subject. 

First, we note that a conductivity formula with a ~ 
— ln 2 (-\/n/no) cannot, by definition, qualitatively ac- 
count for the most important aspect of graphene trans- 
port, namely, the existence of the low-density minimum 
conductivity for a finite range of density around the Dirac 



(i.e. charge neutrality) point. Thus, the resonant scatter- 
ing model, even in the most favorable circumstances, can 
only be a rather phenomenological data-fitting scheme for 
<x(n) in an intermediate density range n c < n < hq where 
n c defines the density regime for the graphene minimum 
conductivity plateau around the Dirac point (taken to 
be at n = 0) and n = (nR 2 )' 1 . Second, as no phys- 
ical evidence for the existence of these "resonant scat- 
tering" short-range atomic defects (with ~ 3 x 10 11 
cm" 2 ) is presented in Ref. (ic| except for providing an 
intermediate-density phenomenological fit to their con- 
ductivity data, it is not clear whether the resonant scat- 
tering is the only scattering mechanism for graphene on 
SrTi0 3 . 

One interesting observation in Ref. fiol ] is the apparent 
absence of much temperature dependence in the graphene 
conductivity while the substrate dielectric constant k(T) 
is changing substantially as a function of temperature as 
the substrate material, i.e. SrTiOa undergoes a para- 
electric to ferroelectric transition with the lowering of 
temperature. One possibility that cannot be ruled out in 
this context is that the substrate impurity density is also 
changing substantially during this temperature-induced 
substrate structural ferroelectric transition. In the cur- 
rent work, we explicitly incorporate such a possibility 
in the theory for disorder limited graphene conductiv- 
ity to compare with the data of Ref. [13]. We believe 
that the presence of some background Coulomb disorder 
is essential for understanding the low-density minimum 
conductivity behavior near the charge neutrality point in 
graphene on any substrate. 

In this paper, we show our best theoretical fits to the 
data of Ref. [l0| in Fig. Q] finding that a dual model 
involving both Coulomb and resonant scattering disor- 
der can well explain the data of Ref. [Io| with the single 
assumption of a variable background charged impurity 
density with varying temperature. The assumption of 
a temperature-dependent charged impurity density for 
SrTiOa substrates is not an arbitrary data fitting ploy 
because the complicated lattice ferroelectric properties 
of SrTiOa leading to the strong functional dependence 
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of the dielectric constant on temperature may very well 
also produce a temperature dependent charged impurity 
density increasing strongly with decreasing temperature 
just as the actual carrier density in graphene on SrTi03 
increases rapidly with decreasing temperature at a fixed 
gate voltage (see Fig. 2 in the Supplementary Informa- 
tion of Ref. 

In the remaining of this paper, we demonstrate that 
resonant scattering by itself can never explain the exper- 
imental data in the whole carrier density range for any 
substrate since the low-density minimum conductivity 
behavior cannot be understood based on a model which 
includes only resonant scattering. We present our best 
theoretical fits to the experimental data in Ref. for 
SrTiOs using the new theory of long-range Coulomb dis- 
order and the resonant scattering defects. We establish 
that a simple model of single type of scattering mech- 
anism can not describe the transport data both at low 
and high carrier density for complex oxide substrates. 
We also compare three different models to explain trans- 
port properties of graphene on different substrates, find- 
ing that different behaviors of cr(n) for different samples 
on different substrates are attributed to the competition 
among different scattering sources. The main message of 
our paper illustrated through concrete calculations, not 
emphasized explicitly in earlier works (although it might 
have been implicit), is that the high-density transport 
in graphene is nonuniversal and reflects the combina- 
tion of various scattering processes arising from the sub- 
strate whereas the low-density transport is always dom- 
inated by long-range Coulomb disorder. In addition to 
the SrTiOs substrates used in Ref. we also consider 
graphene transport on several other common substrates 
(e.g. SiC>2, h-BN, and vacuum) to make our point about 
the universality (nonuniversality) of the transport behav- 
ior at low(high) densities. 
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FIG. 1. Fits to the experimental data of graphene on SrTi03 
(Fig. 2 of Ref. QU) including Coulomb impurity ni mp and 
resonant impurity rii . (a) Resistance p as a function of carrier 
density n. (b) Conductance it as a function of n (in semi- 
logarithmic scale), (c) Conductance a as a function of n (in 
linear scale). The solid (dashed) line is for the temperature 
T = 0.25 (50) K, k = 3500 (500) and n imp = 3 x 10 14 (2 x 10 13 
cm -2 ). Note that we use the theory presented in Ref. Q, s 
denotes the potential fluctuation associated with the puddles 
induced by Coulomb disordei^, the average distance of the 
charged impurity from the graphene sheet used in (a), (b) and 
(c) is d = 1 A and the range of the resonant scattering defects 
R = 0.22 nm. Solid (red) dots represent the experimental 
data extracted from Ref. [IXj. 



conductivity is strictly linear in carrier density for scat- 
tering by long-range Coulomb disorder. 



II. THEORETICAL FORMALISM 

In this section, we first discuss the effects of vari- 
ous scattering mechanisms on the electronic transport of 
graphene. We will analyze four different sources of disor- 
der covering the various realistic possibilities for resistive 
disorder scattering in graphene (a) randomly distributed 
long-range screened Coulomb impurity^, (b) short-range 
point defects^, (c) resonant disorder—, and (d) correlated 
long-range Coulomb impurity^ in the following. The in- 
verse of the total scattering time l/T tot is the sum of 
different inverse scattering times due to various scatter- 
ing mechanisms at low temperatures. We further discuss 
three different models leading to the sub-linear behav- 
ior of cr{n) of graphene on various substrates: (1) the 
standard model including randomly distributed Coulomb 
impurity and short-range scattering mechanism; (2) the 
resonant scattering mechanism; (3) the recently proposed 
correlated charged impurity models. We note that at 
densities above the puddle-dominated density regime, the 



First, the scattering time due to randomly distributed 
charged impurity, denoted as Ti mp , is given by£, 
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x [1 - cos6> kk /] 5(e p v - e pk ) (1) 



where e pk = phvpk is the energy of non- interacting Dirac 
fcrmion with the Fermi velocity vf for the pseudospin 
state "p" and 2D wave vector k, # kk / is the scattering 
angle between in- and out- wave vectors k and k', q = 
k — k', g(# kk ') = [1 + cos 6* kk <] /2 is a wave function form 
factor associated with the chiral matrix of monolayer 
graphene (and is determined by its band dispersion rela- 
tion) . iii mp is the 2D density of the randomly distributed 
screened Coulomb impurity^. Vi(q) = 2ive 2 / (nq) is the 
Fourier transform of the two-dimensional (2D) Coulomb 
potential in an effective background dielectric constant 
k, e(q) is the static dielectric function of graphene within 
random-phase approximation (RPA)^£. This long-range 
Coulomb disorder leads to the linear density dependent 
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conductivity at low temperatures^. 



h 2n imp r 2 s Gi(r s ) 



(2) 



where r s — e 2 / (Kvfk) is graphene fine structure constant 
andGi(a;) = j + 6x — 6ttx 2 + 4x(6x 2 — l)g(x) with g (x) = 
scch _1 (2a-)/Vl - 4a; 2 for x < \ and sec _1 (2a;)/V4x 2 - 1 
for x > \. 

We then provide the short-range disorder scattering 
time T sd , given by& 



k 



n sd V 



(3) 



where n sc i is the 2D short-range impurity density and Vq 
is a constant short-range (i.e. a 5-function in real space) 
potential strength. The conductivity at low temperature 
induced by this short-range impurity has the following 
form, 
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FIG. 2. (Color online). Theoretical fits to a(n) for graphene 
on other substrates. The solid line is the fit including 
Coulomb impurity n imp and short-range point defect n B dVo . 
The dashed line is the fit including Coulomb impurity m m p 
and resonant impurity rn. (a) Fits to the o"(n) data for 
graphene on Si02 (Fig. 2 of Ref. Q]). The potential fluctu- 
ation s = 55 meV, the effective dielectric constant n — 2.45, 
and the average distance of the charged impurity from the 
graphene sheet is d = 1 A. (b) Fits to the er(n) data of sus- 
pended graphene (Fig. 3(c) of Ref. Q3]). The potential fluc- 
tuation s = 10 meV, the temperature T = 5 K, the dielectric 
constant k = 1.0, and the average distance of the charged 
impurity from the graphene sheet is d = 1 A. 



which is independent of carrier density (cr(n) ~ con- 
stant). 

The scattering time due to resonant defects has been 
shown to have the following forroii 
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where rii is the concentration of the resonant defects and 
R is the potential range of the resonant scattering de- 
fects. The logarithmic term in the scattering time of 
resonant defects gives rise to the sublinear density de- 
pendent conductivity 2 ^, 



2e 2 n 
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We now turn to the discussion of the scattering time 
due to the spatially correlated charged impurity r c is 
given by-H 
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where 5(g) = 1 — 27rn im „ — Ji(qr ) is the structure factor 

q 

using the simple continuum analytic models, J\ (x) is the 
Bessel function of the first kind and the correlation effects 
is defined by the length scale r < r t = (7rn imp )~ 1//2 , 
the so-called correlation length. The asymptotic form of 
conductivity at low u kp " arising from correlated long- 
range Coulomb disorder is found to beH, 
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h 2ni mp r 2 Gi(r s ) 1 — a + Ba 2 n/n 
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Tm imp r 2 and B = G 2 (r s )/[2Gi(r s )] with 
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where a 

G 2 (y) = 

correlated Coulomb disorder also leads to a sub-linear 
density-dependent conductivity, which was adopted to 
explain the enhancement of both sublinearity of (j{n) and 
mobility reported in Ref. [i~6j . 

Comparing different scattering mechanisms discussed 
above, it is obvious that the long-range Coulomb disor- 
der dominates <r(n) at low carrier density since the corre- 
sponding scattering rate is asymptotically the largest in 
the vanishing density limit, implying that as the density 
decreases, only Coulomb disorder would matter near the 
Dirac point. The long-range Coulomb disorder is also 
essential in order to explain the electron-hole puddle in- 
duced low-density minimum conductivity^. In addition 
to the resonant and correlated long-range Coulomb scat- 
tering, we note that a standard model involving both 
long-range (a(n) ~ n) and short-range scatterers can 
explain the sublinearity of cr(n) at higher carrier den- 
sity, which has successfully explain the conductivity of 
graphene on boron nitride substrates^. We argue that 
the sample dependent <r(n) at high carrier density may 
well arise from the competition between different scatter- 
ing mechanisms. 



III. NUMERICAL RESULTS 

In this section, we show our theoretical fits to some rep- 
resentative u(n) data for graphene on four different sub- 
strates. The analytical formula for a(n) shown in the last 
section are calculated for homogeneous systems. To cap- 
ture the inhomogencous nature of the graphene landscape 
close to the charge neutrality point (CNP), we assume a 



4 



Gaussian form of potential fluctuation, parametrized by 
"s" and apply the effective medium theory of conduc- 
tance with binary mixture of component, i.e., electron- 
hole puddles, the details of which have been given in 
Ref. [7|. The minimum conductivity plateau of the pud- 
dle dominated regime in graphene around the charge neu- 
trality point is just well-captured in our theory. 

In Fig. [TJ we show our best theoretical fits to graphene 
transport data of Ref. [l(| , using a dual scattering mod- 
els: the new theory of long-range Coulomb disorder 
and the resonant scattering defects, which give cr(n) ~ 
n\\^{\rmB?). The Coulomb disorder theory of trans- 
port follows Ref. and includes the effect of Coulomb 
disorder-induced electron-hole puddles through the po- 
tential fluctuation parameter "s" . The new model ex- 
plains the conductivity of graphene on SrTi03 in the 
whole range of carrier density. Couto et al. in Ref. 
pointed out that at fixed value of external gate voltage, 
the carrier density n and the dielectric constant of SrTi03 
substrate would increase by approximately one order of 
magnitude as the temperature is lowered from 50 K down 
to 250 mK, which may also lead to the temperature- 
dependent charged impurity density. Resonant scattering 
may account for the conductivity of graphene on SrTi03 
away from the Dirac point. Note, however, that Coulomb 
disorder induced electron-hole puddles, which give rise 
to the existence of minimum conductivity, dominate the 
transport properties of graphene on SrTi03 close to the 
Dirac point and thus some charged impurity scattering 
must be present in the samples of Ref. ■ 

In Fig. [21 we present the comparison of the er(n) re- 
sults between the standard model and the new model 
for graphene on Si02 substrate 2 - and the suspended 
graphen o 17 ' 18 . Both models can fit the experimental data 
very well. Fig.[5]shows that the above mentioned two dis- 
tinct dual scattering models are equally successful in de- 
scribing the experimental data for both Si02-bascd and 
suspended graphene, and there is no particular reason, 
within the transport data fitting scheme without addi- 
tional information, to preferably choose one model over 
the other. We mention that the parameters of the models 
(n-imp, n s dV^ in the standard model and Tiimp j f^i in the 
new model) are not absolutely unique, and it is possi- 
ble to get equivalent fits by adjusting the parameter sets 
somewhat. This is expected because cr(n) is a smooth 
function at higher density and somewhat different pa- 
rameter sets for disorder cannot be distinguished since 
the models have only qualitative and semi-quantitative 
predictive power. We emphasize, however, that the res- 
onant scattering model by itself cannot explain the min- 
imum conductivity phenomenon around the charge neu- 
trality point in graphene, existing all the way to the room 
temperature and above, which is a generic observation for 
graphene on all substrates. 

In Fig. [3ja) , we compare the theoretical fits between 
the new model and the correlated disorder model pro- 
posed in Ref. [l3| for suspended graphene. Because 
the annealing procedure is routinely used in preparing 
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FIG. 3. (Color online), (a) Fits to the <r(n) data of suspended 
graphene (Fig. 3(c) of Ref. [TJ]). The solid line is the fit 
including the Coulomb impurity = 0.15 x 10 cm and 
resonant scattering m. The dashed line is the fit including 
correlated Coulomb impurity ni mp = 1.6 x 10 10 cm -2 and 
correlation length ro = 60ao with ao = 4.92 A. (b) Fits to 
the crjn) data for graphene on boron nitride (Fig. 3(a) of 
Ref. Q). The solid line is the fit including Coulomb impurity 
Tiimp = 1.4 x 10 11 cm -2 and short-range point defect n a dV§ , 
the potential fluctuation s = 44 meV. The dashed line is 
the fit including correlated Coulomb impurity rii mv = 3.7 x 
10 11 cm~ 2 , and the potential fluctuation s = 38 meV. Note 
that the average distance of the charged impurity from the 
graphene sheet used in (a) and (b) is d = A. 



suspended graphene samples, it is very likely to intro- 
duce some spatial correlation among the Coulomb impu- 
rities. In Fig. [3Jb) , we compare the standard model with 
the correlated disorder model for graphene on hexagonal 
boron nitride (h-BN) substrate. Graphene on boron ni- 
tride may also induce some inter-impurity correlations 
imposed by the similarity between boron nitride and 
graphene lattice structure. We know, from Fig. [3J all of 
these three models can give reasonable fits to the experi- 
mental data of the density-dependent conductivity. They 
capture the experimentally observed strong-nonlincarity 
in the density-dependent conductivity a(n). The long- 
range Coulomb impurity leads to the density inhomo- 
geneity in the graphene system, which is essential in ex- 
plaining the conductivity plateau close to the Dirac point. 



IV. DISCUSSION AND CONCLUSION 

Before concluding, we argue that the resonant scat- 
tering theory with a ~ — In 2 ( y/n/no) predicts a van- 
ishing of graphene conductivity as the carrier density 
approaches the Dirac point (n = 0) which is not ob- 
served in the experiments of Ref. [lOj, where cr(n) be- 
comes an approximate constant at low density around 
the Dirac point in complete contrast with the resonant 
scattering predictions and thus, the resonant scattering 
theory cannot obviously be the complete story underly- 
ing the transport mechanisms in Ref. ■ We emphasize 
that the high-density behavior of er(n) in graphene is, in 
general, knowr£ to be nonuniversal with different sam- 
ples on different substrates showing different behaviors^ 
arising from the competition among Coulomb disorder, 
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short-range disorder, ripples, and (perhaps even) reso- 
nant scattering disorder. It is entirely possible that the 
transport data of Ref. is best described by a com- 
bination of Coulomb disorder and resonant scattering, 
where the low-density minimum conductivity arises from 
the Coulomb disorder and the intermediate (sublinear 
in) density conductivity arises from resonant scattering 
as shown in FigfTJ assuming dual independent scatter- 
ing by charged impurity and resonant scattering centers. 
The fact that we get excellent agreement over the whole 
density range of the experimental data indicates that 
our dual scattering model is a more reasonable descrip- 
tion than the pure resonant scattering model. We have 
also obtained similar good theoretical fits to the existing 
graphene data as shown in Fig. [2] and [3] Thus, the pos- 
sibility that the resonant scattering mechanism is opera- 
tional at some level in graphene transport at higher den- 
sities on any substrate (not just SrTiOa) cannot defini- 
tively be ruled out as we show in this paper by consid- 
ering several different substrates. What can be stated 
rather definitively is that the low-density transport must 
always be dominated in graphene by charged impurity 
scattering. 

We conclude by emphasizing that the resonant scatter- 
ing model by itself can not explain the transport data in 



graphene. In particular, the long-range charged impurity 
scattering is needed in order to explain the existence of 
the minimum conductivity close to the charge neutral- 
ity point. The main point of our paper is that the un- 
derstanding of the graphene minimum conductivity phe- 
nomenon necessarily requires the inclusion of scattering 
by random charged impurities in the environment. We 
emphasize that since a r /<7 c ~ ln 2 (-^/ njno), where a r<c 
are the resonant scattering and Coulomb scattering in- 
duced conductivity respectively, Coulomb disorder—, with 
<7 C ~ n , must necessarily dominate graphene resistivity 
as one approaches the n — > charge neutrality point, 
not just in SrTiC>3, but in all systems. The high density 
graphene transport is in general, however, nonuniversal 
where many scattering mechanisms may contribute de- 
pending on the density range. 
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